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1. Introduction 

For the description of nuclear processes in many astrophysical scenarios the knowledge of 

a-nucleus potentials is necessary. Such processes are radiative capture, transfer reactions 

Q\ \ and alpha-decay occuring in primordial nucleosynthesis, stellar hydrostatic and explosive 

burning modes. 
5h ; 

Until now such nuclear processes have often been described using mainly phenomenological 
and energy-independent potentials (e.g. square well, Woods-Saxon potentials etc.). In this 
work we develop a-nucleus potentials using the folding procedure [1,2]. With this method 
the ambiguity of the phenomenological potentials can be avoided to a great extent. The 
uniqueness and the energy dependence of these potentials are an important feature with 
respect to astrophysical applications. 

(N : 

Such a-nucleus potentials have been used successfully for the description of scattering 
processes [3,4,5], transfer reactions [6,7,8] and radiative capture [9,10] on light nuclei. In 
this work we extend our systematic investigation of a-nucleus potentials to intermediate 
and heavy stable and unstable nuclei (A > 70) which are relevant for the p-process [11]. 
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2. Folding procedure 

The real part of the optical potential is deduced in the framework of the double-folding 
model of Kobos et al. [1] and is described by 

V(r) = X J d?i J dr 2 pT(ri)Pa(r2)t(E,p T ,p a ,s = r + f2-fi) , (1) 

where r is the separation of the centers of mass of the colliding target nucleus and the a 
particle, prifi) and p a (r2) are the respective nucleon densities derived from nuclear charge 
distributions [12], and t(E, pt, p a , s) is the density-dependent effective NN interaction [1]. 
By means of the normalization factor A the depth of the potential can be adjusted to 
elastic scattering data and to bound and resonant state energies of nuclear cluster states. 

The imaginary part of the potential can be parametrized by a Woods-Saxon form. 

The strength of the potential is measured by its volume integral per interacting pair of 
nucleons, e.g. for the real part 

J R {E) = — \ V(r,E)r 2 dr , (2) 
Ap ■ At Jo 
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Fig. la: Volume integrals of the real part 
of the optical potential for the a + 4 He, 
a + 16 O, and a + 40 Ca system. 



Fig. lb: Volume integrals of the imagi- 
nary part of the optical potential for the 
a + 4 He, a + 16 0, and a + 40 Ca 



system. 

where Ap and At denote the projectile and target mass numbers, respectively. 
In Figs, la and lb the volume integrals for the real part, Jr, and the imaginary part, 
J/, of the optical a-nucleus potential for the target nuclei 4 He, 16 and 40 Ca are shown. 
These values have been obtained from the analysis of elastic scattering data and from 
calculations of 8 Be, 20 Ne and 44 Ti cluster states, respectively [3,4]. A strong energy and 
mass dependence for both, the real and imaginary parts, can be observed. The energy 
dependence of Jr (Fig. la) is due to the energy dependence of both, the effective NN 
interaction (Eq. 1) and the so-called dynamic polarization potential, which is related to 
that of the imaginary part of the potential by a dispersion relation [13,3]. The curves 
shown in Fig. la represent the results of calculations which consider both effects [4]. In 
these calculations the energy dependence of Jj was parametrized by [14] 



with Eq being the threshold energy for inelastic processes. A linear regression procedure to 
the data points results in values for Jq and A. The curves calculated with these parameters 
are given in Fig. lb as solid, dashed and dotted lines. 

3. Alpha-nucleus potentials for intermediate and heavy nuclei 

Using the same procedure as described in Sec. 2, we determined optical alpha- nucleus 
potentials for intermediate and heavy nuclei. In Fig. 2 calculated elastic scattering cross 
sections for some target nuclei are compared with experimental data at different energies. 
An excellent agreement between the theoretical analysis and the experiment is obtained. 
From this fit to the experimental data, the strengths of the real part, A, as well as the 
parameters of the imaginary Woods-Saxon potentials are obtained for energies above 15 






MeV. 
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Fig. 2: Elastic a scattering on Zr, Pb, Ge, 141 Pr and Sm: Experimental data 
[25-37] and optical model fits calculated by using double-folded potentials. 
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Fig. 3a: Volume integrals of the real part 
of the optical ct-nucleus potential for some 
intermediate and heavy nuclei. 



Fig. 3b: Volume integrals of the imag- 
inary part of the optical a-nucleus po- 
tential for some intermediate and heavy 
nuclei. 



The volume integrals for the real as well as the imaginary part of these potentials are 
shown in Figs. 3a and 3b. The following results can be deduced: 

(i) The mass dependence of the volume integrals for the real part of the potentials is very 
weak for all heavier nuclei. However, the absolute values are somewhat smaller than 
those for lighter nuclei (compare Figs, la and 3a). 

(ii) The energy dependence has a similar form as the one obtained for lighter nuclei. The 
volume integrals of the real part have a maximum about Ecm = 30 MeV and decrease 
slightly when going to lower and higher energies. For astrophysically relevant energies, 
which are in the order of 10 MeV, a linear extrapolation to lower energies is performed 
(s. Fig. 3a). We obtain 



The curve shown in Fig. 3a is the result of a calculation for the a- Pb potential 
which contains the energy dependence of the folding potential and the dispersive 
part. The knowledge of bound and quasibound-state potentials is also necessary for 
the calculation of transfer and capture cross sections. We have calculated these alpha- 
nucleus potentials for some bound states. The volume integrals for these potentials 
are also close to the extrapolated values (s. Fig. 3a). 

(iii) The volume integrals for the imaginary part of the potentials obtained from the fit 
to the experimental elastic scattering data for the target nuclei 70 Ge, 90 Zr, 141 Pr, 
144 Sm and 208 Pb are shown in Fig. 3b. For 208 Pb the parametrization given in Eq. (3) 
was used to calculate the observed energy dependence of Jj (solid line in Fig. 3b). 
As expected a strong mass dependence of the volume integrals is observed, since the 
strength of the imaginary potential is quite different for a doubly magic and a strongly 
deformed nucleus. Therefore for an energy of Ecm = 10 MeV, J/ values can range 
between 10 and 60 MeV fm 3 . 



Jr/(A p ■ A T ) = (320 + 0.67£ CM ) MeV fm 3 

(E C M in MeV) 



(4) 
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For unstable nuclei the mass densities necessary for the calculation of the folding potential 
cannot be obtained from electron scattering data. In these cases the densities can be 
calculated in the (crwp) model [15,16] used in the relativistic mean-field theory. We used 
the parameter set NLSH which is suited for neutron and proton rich nuclei [17,18,19]. We 
found that for stable tin and samarium isotopes the densities calculated in this model 
compare well with the experimental data [12]. In order to calculate the alpha-nucleus 
potentials for unstable nuclei, the strengths of the potentials were adjusted to reproduce 
the parametrized volume integrals given above. 

4. Application to "p-process" isotopes 

The alpha-nucleus potentials determined with the folding procedure are necessary for the 
calculation of (7, a) photo disintegration cross sections in the p-process. As an example 
we consider the inverse reaction 144 Sm(a, 7) 148 Gd. This reaction determines the ratio 
142 Nd/ 144 Nd in some meteorites [20]. In previous work the astrophysical S-factors and 
photonuclear reaction rates have been generated using the statistical theory of nuclear 
reactions as employed by Michaud and Fowler [21]. They used an equivalent square well 
[ESW] deduced from a Woods-Saxon potential. However, the effective radius parameter 
for this ESW is quite uncertain [22]. For the two different ESW-radii for the a- 144 Sm 
potential (R a = 8.75 fm and 8.01 fm), the calculated cross sections differ by a factor of 
ten. 

Table 1: S-factors and reaction rates for 144 Sm(«, 7) 148 Gd 





ESW 
R = 8.75 fm 


ESW 
R = 8.01 fm 


Woods-Saxon 
potential [24] 


Folding 
potential") 


S- factor 

(E a = 9.5 MeV) [MeV • b] 


2.3 • 10 28 ») 


2.3 • 10 27 


1.2 • 10 28 


7.8 • 10 27 


reaction rate 

T 9 = 2.5 [cm 3 mol" 1 s" 1 ] 


3.75 • 10" 15 


3.72 • lO" 16 


1.95 • 10" 15 


1.27- 10" 15 


reaction rate 

T 9 = 3.0 [cm 3 mol" 1 s" 1 ] 


2.35 • 10" 12 


2.58- 10" 13 


1.22 • 10" 12 


7.56- 10" 13 



a ) this work 

b ) There is a misprint in [22] giving this value as 2.3 ■ 10 



In table 1 we list the astrophysical S-factor at the Gamow energy E a = 9.5 MeV and the 
reaction rates T 9 = 2.5 and T 9 = 3.0. In the first two columns the results of calculations 
using ESW [22] are given. In the third and fourth column the results of Hauser-Feshbach 
calculations using the code SMOKER [23] are shown. In the first case (column 3) an 
energy-independent Woods-Saxon potential with V = 185 MeV, W = 25 MeV, R = 
R w = 1.4- A 1/3 fm and a = aw = 0.52 fm [24], in the second case (column 4) for the real 
part a folding potential (A = 1.1573) and for the imaginary part a Woods-Saxon potential 
(W = 10 MeV, R w = 1.4 • A 1/3 fm, a w = 0.52 fm) was used. 
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With our improved folding potential the reaction rates shown in the last column of table 
1 are about 1/3 of the value for an ESW with a radius of 8.75 fm. This corresponds about 
to 1/3 of the "recommended" value for the reaction rate in [22] giving a similar reaction 
rate as case C in table 1 of [22]. Therefore, our 146 Sm/ 144 Sm ratio is about 0.22 which is 
consistent with the cosmochemical data of 0.1 - 0.7 [20]. 

The astrophysical S-factor at E a = 9.5 MeV calculated in the direct-capture model gives 
an upper limit of about 10 22 MeV • b, which is more than five orders of magnitude smaller 
than the Hauser-Feshbach result. 
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